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Discoveries of 3 K-shell Lines of Iron and a Coherent Pulsation of 

593-sec from SAX J1748. 2-2808 
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Abstract 

SAX J1748.2— 2808 is a unique X-ray object with a flat spectrum and strong emission lines at 6.4- 
7.0 keV. The Suzaku satellite resolved the emission lines into 3 K-shell lines from neutral and highly ionized 
irons. A clear coherent pulsation with a period of 593-sec was found from the Suzaku and XMM-Newton 
archives. These facts favor that SAX J1748.2— 2808 is an intermediate polar, a subclass of magnetized 
white dwarf binary (cataclysmic variable: CV). This paper reports on details of the findings and discusses 
the origin of this source. 
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1. Introduction 

SAX J1748.2-2808 was discovered with the Beppo- 
SAX satellite in the direction of the radio complex 
Sagittarius (Sgr) D (Sidoli et al. 2001), and was deeply 
re-observed with the XMM-Newton satellite (Sidoli et al. 
2006) . The X-ray spectrum was well- fitted with a power- 
law of photon index 1.4 plus a broad line of 0.43 keV 
(ler) with the center energy at 6.6 keV. The Ah value 
was as large as 1.4 x 10 23 H cm" 2 (Sidoli et al. 2006). 
SAX J1748. 2-2808 was thus regarded as one of the bright- 
est samples of resolved point sources in the Galactic cen- 
ter region (GC) with Chandra deep-exposure observations 
(Muno et al. 2003, Muno et al. 2006). The integrated 
spectra of the point sources resemble the Galactic cen- 
ter diffuse X-rays (GCDX) in the iron line features close 
to 6.4-7.0 keV (Muno et al. 2004, Koyama et al. 2007b). 
Therefore, these, including SAX J1748. 2-2808, can be 
regarded as significantly contributing to the GCDX, even 
though SAX J 1748. 2 -2808 was brighter than most of the 
other point sources in the GC. 

The XMM-Newton spectrum favored that SAX 
J1748.2-2808 is a High Mass X-ray Binary (HMXB) pul- 
sar located at the GC region (Sidoli et al. 2006). However, 
no coherent pulsation was reported from either the Beppo- 
SAX or the XMM-Newton observations. Furthermore, the 
broad line of 0.43 keV (ler) is very unusual as a HMXB. If 
this broadening is due to the Doppler effect, the velocity 
dispersion is as large as 2 x 10 4 km s _1 . No such HMXB 
has been reported so far (e.g. Nagase 1989). The detailed 
spectrum of SAX J1748. 2-2808 did not indicate whether 
the spectral nature is thermal or non-thermal. 

In order to understand the nature of this unique ob- 
ject, we analyzed the Suzaku data of two pointing obser- 
vations on the supernova remnants (SNR), Sgr D SNR 
and G 0.9+0.1, and also re-analyzed the XMM-Newton 
archive data. 



2. Observations and Data Reduction 

The Suzaku and XMM-Newton observations of SAX 
J1748. 2-2808 arc listed in tabic 1. The Suzaku satel- 
lite observed the two SNRs, Sgr D SNR and G 0.9+0.1, 
where SAX J1748. 2-2808 was located near the edge of 
the field of view (FOV) of the XIS. The XIS consists of 
four sets of X-ray CCD camera systems (XIS 0, 1, 2, and 
3) placed on the focal planes of four X-Ray Telescopes 
(XRT) aboard the Suzaku satellite. XIS 0, 2, and 3 
have front-illuminated (FI) CCDs, while XIS 1 has a 
back-illuminated (BI) CCD. One of the FI CCD cam- 
eras (XIS 2) has been out of function since September, 
2006, and hence we did not use it. Detailed descriptions 
of the Suzaku satellite, the XRT, and the XIS can be 
found in Mitsuda et al. (2007), Serlemitsos et al. (2007), 
and Koyama et al. (2007a). 

The XIS observations were made with the normal mode. 
The XIS pulse-height data for each X-ray event were con- 
verted to Pulse Invariant (PI) channels using the xispi 
software version 2007-03, and the calibration database 
version 2008-01-31. We removed the data during the 
epoch of low-Earth elevation angles less than 5 degrees 
(ELV< 5°), day Earth elevation angles less than 10 degrees 
(DYE_ELV< 10°), and the South Atlantic Anomaly. The 
good exposure times are listed in table 1. Although the 
XIS CCDs were significantly degraded by on-orbit par- 
ticle radiation, the CCD performances were restored by 
the Spaced-low charge injection technique (Uhicyama et 
al. 2009). Then, the overall spectral resolutions at 5.9 
keV were -150 and ~170 eV (FWHM) for the FI and 
BI CCDs, respectively. We analyzed the data using the 
software package HEASoft 6.4.1. For the spectral fittings, 
we made XIS response files using xisrmfgen, and auxil- 
iary files using xissimarf gen. Since the spectrum of the 
non- X-ray background (NXB) depends on the geomag- 
netic cut-off rigidity (COR) (Tawa et al. 2008), we ob- 
tained COR-sorted NXB spectra using xisnxbgen, from 
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Table 1. Observation data list 
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Obs. ID 


Date 
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fvvvv / mm / dd 1 


(ksec) 


Suzaku/XIS 


Sgr D SNR 


502020010 


2007/09/06 


139.1 


Suzaku/XIS 


G 0.9+0.1 


502051010 


2008/03/11 


138.8 


XMM-Ncwton/MOS 


Sgr D SNR 


0112970101 


2000/09/23 


15.7 


XMM-Newton/PN 


Sgr D SNR 


0112970101 


2000/09/23 


11.5 


XMM-Newton/MOS 


G 0.9+0.1 


0144220101 


2003/03/12 


49.5 


XMM-Newton/MOS 


SAX J 1748. 2 -2808 


0205240101 


2005/02/26 


50.2 


XMM-Newton/PN 


SAX J1748. 2-2808 


0205240101 


2005/02/26 


41.5 



After the data screening described in the text. 



the night-Earth data released by the Suzaku XIS team. 

SAX J1748.2-2808 was also in the filed of view of the 
XMM-Newton observations on Sgr D SNR and G 0.9+0.1 
on September 2000 and March 2003, respectively. The 
pointing observation on SAX J1748.2— 2808 was also made 
on February, 2005. The X-ray data were obtained with 
the European Photon Imaging Camera (EPIC) (Striider 
et al. 2001; Turner et al. 2001) in an extended full-frame 
mode. The data were analyzed using the Science Analysis 
Software (SAS 7.1.0). Event files for both the PN and the 
Metal Oxide Semiconductor (MOS) detectors were pro- 
duced using the epchain and emchain tasks of SAS, re- 
spectively. The event files were screened for high particle- 
background periods. Good exposure times arc listed in 
table 1. In our analysis, we used the events corresponding 
to patterns 0-4 for the PN and 0-12 for the MOS instru- 
ments. The PN data of SAX J1748. 2-2808 on the 2003 
observation was not available. 

3. Analysis and Results 

3.1. X-ray Image 

The X-ray images taken with the XIS were analyzed 
for each observation separately. Since the Suzaku nom- 
inal position error is < 20" (Uchiyama et al. 2008), we 
fine-tuned the Suzaku position while referring to the po- 
sitions of XMM-Newton sources. In the Sgr D SNR ob- 
servation, we found an XMM-Newton source S10 (see ta- 
ble 4 of Sidoli et al. 2006). The XMM-Newton position 
of S10 is (a, 8) = (267.°24531, -28.°188918), while the 
Suzaku position is (a, 6) = (267.°24409, -28.°185869). 
We therefore shifted the Suzaku coordinate by (Aa, AS) 
= (0.°00122, -0.°003049). In the GO.9+0.1 observa- 
tion, we used an XMM source S2 (see table 4 of 
Sidoli et al. 2006). The XMM Newton position is 
(a, 5) = (266.°81562, -28.°181511), while the Suzaku 
position is (a, 5) = (266.°81737, -28.°180065). We 
hence shifted the Suzaku coordinate by (Aa, A<5) = 
(-0.°00175, -0.°001446). After these fine-tunings of the 
Suzaku coordinate, we made the X-ray image shown in fig- 
ure 1. We found two sources near to the edge of each field. 
These two sources coincide in position with S3, i.e., SAX 
J1748.2— 2808 (stronger source) and S12 (fainter source) 
of the XMM-Newton observation within statistical errors 



In order to estimate the intensity ratio of the two 
sources, we made a projected profile along the line con- 
necting SAX J1748.2-2808 and S12 (see figure lc). The 
profile was fitted with two Gaussians plus a liner function, 
as is shown in figure 2. The widths (1-cr) of the Gaussians 
approximate the projected point-spread function. The 
best-fit source fluxes (normalizations of the Gaussians) 
were determined to be 7.5 x 10 -3 and 2.8 x 10 -3 counts 
s" 1 , for SAX J1748. 2-2808 and for S12, respectively. 

3. 2. Spectrum 

Since the X-ray spectrum of XMM-Newton has al- 
ready been reported, we show the Suzaku spectrum of 
SAX J1748. 2-2808. The point-spread function of the 
XRT has a complicated polygon shape images, due to- 
gether to the 4-segments structure of the XRT and image 
deformation near the field edge (Serlemitsos et al. 2007). 
We therefore extracted the spectra from the solid poly- 
gons in figures la and lb, for the best S/N ratio within 
limited statistics. 

The background spectra were obtained from a near-by 
sky shown by the dashed squares, from which the dashed 
elliptical regions were excluded. The background spectra 
consist of the non- X-ray background (NXB) and the lo- 
cal background (the cosmic X-ray background plus the 
Galactic center diffuse X-rays). Since the local back- 
ground is affected by the vignetting of the XRT while the 
NXB is not, each background was treated separately. For 
both the source and the background data, we made COR 
distributions and composed the NXB spectra from the 
COR-sortcd NXB data set (see section 2). We then sub- 
tracted the NXB extracted from the same detector areas. 
After NXB subtraction, we corrected the vignetting effect 
due to the different off-axis angles between the source and 
the background regions by multiplying the effective-area 
ratios for each energy bin of the local background spec- 
tra (the same process as Hyodo et al. 2008). We further 
subtracted the corrected local background from the source 
spectra. All of the source spectra of the two FIs from the 
two observations are co-added to increase the statistics. 

We estimated the contamination of the near-by fainter 
source S12 by a ray-tracing simulation (xissim in 
the HEASoft package). The spectral shape of SAX 
J1748.2— 2808 was assumed to be an absorbed power-law 
plus 3 Gaussian lines, whose parameters are the same as 
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Fig. 1. X-ray image in the 3—7 keV band from the Suzaku observations on Sgr D SNR (a) and G 0.9+0.1 (b). The source regions 
of SAX J1748.2— 2808 are given by the solid polygons, which trace the complicated point-spread function of the XRT in the field 
edge (see text). The dashed squares (excluding ellipses) are the local backgrounds, (c): Combined images of figures (a) and (b) near 
SAX J1748.2— 2808. The projected profile (figure 2) was made from the dashed rectangle region, where the white circles indicate 
the positions of SAX J1748. 2-2808 (stronger source) and S12 (fainter source) (Sidoli et al. 2006). 



table 2. This spectral model will be mentioned in the 
next paragraph. Since the fainter source S12 is too weak 
to allow a spectral analysis, its spectral shape was as- 
sumed to be the same model of an absorbed black-body 
(N H = 10 24 cm" 2 , kT =1 keV) as Sidoli et al. (2006). 
Using the flux ratio of SAX J1748. 2-2808 and S12 de- 
termined in section 3.1, we found that the contamination 
was ^$10% in the 3-10 keV band. We therefore ignore the 
contamination of S12 in the following analysis and discus- 
sion. 

The X-ray spectra of FI and BI were simultaneously 
fitted with an absorbed power-law plus one broad line, 
Absx (Power-Law + 1 Gaussian line), which is the same 
model as that of the XMM-Newton by Sidoli et al. (2006). 
The best-fit (90% confidence range) photon index of the 
power-law, line center energy, and width (la) are 1.2 (0.8- 
1.5), 6.66 (6.54-6.77) keV, and 0.31 (0.21-0.39) keV, re- 
spectively (x 2 /dof of 72.8/64). These parameter values 
are consistent with those of XMM-Newton within the sta- 
tistical error. However, the Suzaku spectrum exhibits sig- 
nificant residuals near to energies of 6.4, 6.7, and 7.0 keV. 
We therefore fitted with a model of an absorbed power- 
law plus 3 narrow lines near to 6.40, 6.68 and 6.97 keV, 
Absx (Power-Law + 3 Gaussian lines). Although the fit- 
tings were made simultaneously for the FI and BI spectra, 
we simply show the FI result in figure 3. The fit was im- 
proved with x 2 /dof of 67.6/64. The best-fit parameters 
are listed in table 2. 

The best-fit line energies of 6.7 keV and 7.0 keV are 
likely due to Ka lines from He-like and H-like irons at en- 
ergies of 6.68keV and 6.97 keV, respectively, while the 6.4 
keV line would be a Ka line from neutral irons. The 6.97 
keV line may be contaminated by a K/3 line (7.05 keV) of 
neutral irons. We therefore applied a more physical model: 
an absorbed thin thermal plasma (APEC) plus the 6.40 
keV and 7.05 keV lines from neutral irons, Absx (APEC 
+ 6.40 keV + 7.05 keV lines). In this model, we fixed 
the flux ratio of the 6.40 keV and 7.05 keV lines to 1: 
0.125 (Kaastra & Mewe 1993). A simultaneous fit for the 
FI and BI spectra with x 2 /dof of 77.4/65 is acceptable at 
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Fig. 2. Projected profile along the line connecting SAX 
J1748.2-2808 and S12. This profile was made from the data 
in the dashed rectangle in figure lc. The solid lines are the 
best-fit liner function and 2-Gaussians. 

the 90% confidence level, as is given in table 2. In figure 4, 
we show only the FI spectrum for simplicity. 

3. 3. Timing 

We examined a long-term (8 years) X-ray flux history 
from the 5 observations listed in table 1. The X-ray fluxes 
were calculated by fitting the 3-line model with the param- 
eters fixed values given in table 2. Only the normalization 
(flux) was a free parameter. The resultant fluxes in the 
3—10 keV band show no significant variation during the 8 
years interval. 

We then searched for a coherent pulsation in the 
3—7 keV band from all of the observations listed in table 1. 
The Fast Fourier Transform (FFT) analysis revealed a 
clear peak at ~ 1.7 x 10~ 3 Hz from all of the observations. 
We then searched for accurate pulse period with the fold- 
ing technique, and found a significant peak near to the 
trial period of 593 sec. The best-fit pulse periods and er- 
rors are listed in table 3 for all of the observations. As 
an example, we show the Suzaku results from 2007: the 
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Fig. 3. Background-subtracted spectrum of SAX J1748.2— 2808 for the FI CCD. The source spectrum was made from the solid 
polygons in figures la and lb, while the local background spectra were taken from the dashed squares excluding the dashed ellipses. 
The solid line is the best-fit model of an absorbed power-law plus 3 Gaussians. 
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Fig. 4. Same as figure 3, but the best-fit model of an absorbed thin thermal plasma (APEC) plus K-shell lines of neutral iron. 
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Table 2. Best-fit parameters 
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Model: Abs x (Power-Law + 3 Gaussian Lines) 


Line Energy (keV) Flux (10 b ph s 1 cm 2 ) 


Equivalent Width (eV) 


6.40 (6.39-6.47) 1.5 (0.8-2.3) 
6.68 (6.66-6.72) 1.9 (1.1-2.7) 
6.97 (6.94-7.46) 1.3 (0.6-2.1) 


140 (30-270) 
180 (30-350) 
130 (< 270) 


Parameter Bets-Fit Value 


Photon Index 1.0 (0.8-1.3) 




N R (10 23 cm- 2 ) 1.3 (1.0-1.7) 
Flux (3-10 keV) (10~ 13 erg s" 1 cm" 2 ) 6.0 (5.7-6.2) 
Luminosity* (3-10 keV) (10 33 erg s" 1 ) 7.4 (7.0-7.6) 
X 2 1 dof 67.6 / 64 




Model: Absx(APEC + Neutral Iron Lines) 


Line Energy (keV) Flux (10 b ph s 1 cm 2 ) 


Equivalent Width (eV) 


6.40 (fixed) 2.1 (1.0-2.7) 
7.05 (fixed) 0.26* 


160 (10-270) 
23 


Parameter Best-Fit Value 



Temperature: kT (keV) 
Metal abundance (solar) 



Nn (10 



23 



Flux (3-10 kcV) (10~ 13 erg s" 1 cm- 2 ) 
Luminosity* (3-10 kcV) (10 33 erg s" 1 ) 
X 2 I dof " 



12 (9-17) 
0.57 (0.35-0.85) 
1.8 (1.6-2.1) 
5.4 (5.2-5.7) 
8.4 (8.0-8.8) 
77.4 / 65 



* Parenthesis is 90% confidence range. 

* Absorption corrected. Distance toward SAX J1748.2— 2808 is assumed to be 



* The flux of the 7.05 keV line is constrained to be 0.125 of the 6.40 keV Line. 



.5 kpc. 



power spectrum of FFT, periodogram, and folded light 
curve in figure 5a, 5b, and 5c, respectively. 

The 593-sec pulsation is likely to be a spin rotation of 
either a magnetic white dwarf or a neutron star in a binary 
system (see section 4). We, therefore, searched for an or- 
bital modulation in the light curve with 1000-scc binning. 
In the Suzaku observations, we found no sign of orbital 
modulation, nor any sign of an eclipse. And no orbital 
modulation was found from the XMM-Newton data, ei- 
ther. This is reasonable because the time coverages of the 
XMM-Newton observations were less than the Suzaku. 

4. Discussion 

We found a coherent pulsation of 593-sec from SAX 
J1748.2— 2808. This constrains the origin of this source 
to be either a magnetic CV, or a HMXB pulsar. The 593- 
sec period is among the slowest 10% of the pulse period in 
HMXB pulsars (Nagase et al. 1989, Liu et al. 2000), while 
the fastest 30% in magnetic CVs (Ritter & Kolb 2003). 
Thus, the spin period of 593-sec favors a magnetic CV 
scenario, although a HMXB pulsar scenario is not firmly 
excluded. 

The Suzaku spectrum resolved the broad line at 6.6 
keV found in the previous XMM-Newton observation into, 
at least, three lines at 6.40, 6.68, 6.97 keV. The best-fit 
equivalent width (EW) of these three lines are ~140 eV, 
~180 eV, and —130 eV, respectively. Ezuka and Ishida 
(1999) compiled the ASCA data and reported that the 



magnetic CVs exhibit 3 iron Ka lines at 6.4, 6.7, and 7.0 
kcV with mean EW values of -100, - 200, and -100 eV, 
respectively (see also Hellier et al. 1998), nearly the same 
as those of SAX J1748. 2-2808. 

The overall spectrum was also well fitted with a thin 
thermal plasma of 12-keV temperature with a sub-solar 
iron abundance plus Ka (6.40 keV) and K/3 (7.05 keV) 
lines from neutral irons. Ezuka and Ishida (1999) also 
reported that the spectra of the magnetic CVs can be 
described by a thin thermal plasma model plus 6.4 kcV 
line with a mean temperature of —20 keV, closely resem- 
ble to the Suzaku results of SAX J1748. 2-2808. On the 
other hand, HMXB pulsars show a hard continuum spec- 
trum, but it is a broken power-law and not a thin thermal. 
HMXBs exhibit an iron emission line feature. However, 
the line feature is not complex, but is a single line, mostly 
at 6.4 kcV (see table 3 of Nagase et al. 1989). These 
properties are different from those of SAX J1748.2— 2808. 

All of the above facts favor the idea that SAX 
J 1748. 2 -2808 is a magnetic CV rather than a HMXB. 
The pulse period of 593-sec is smaller than the possible 
orbital period, although we found no orbital modulation. 
Therefore, 593-sec would be the spin period, and hence 
SAX J1748. 2-2808 is an intermediate polar (IP), not a 
polar with synchronized spin and an orbital period. If the 
luminosity of SAX J1748.2— 2808 is typical value for an 
IP of — 10 33 erg s _1 (Patterson 1994), then from the ob- 
served flux of ~ 6 x 10 -13 erg cm -2 s -1 , the distance is 
estimated to be -4 kpc. Thus, SAX J1748.2-2808 would 
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(a) FFT Suzaku (2007) 
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Fig. 5. (a): Power spectrum (FFT) in the 3—7 keV band of the Suzaku observation (2007). (b): Same as (a), but the pcriodogram 
at around the trial period of 593 sec. (c): Same as (a), but the folded light curve at the best-fit period of 593.1 sec, where the dashed 
line is the background level. 



Table 3. Best-fit Period 



Observatory (Year/Month) Instruments Pulse Period (s) 



XMM-Newton (2000/09) 


MOS1- 


h2 


592 


±8 


XMM-Newton (2000/09) 


PN 




594 


±8 


XMM-Newton (2003/03) 


MOS1- 


-2 


593 


±2 


XMM-Newton (2005/02) 


MOS1- 


-2 


595 


±3 


XMM-Newton (2005/02) 


PN 




593 


±2 


Suzaku (2007/09) 


XIS 0+1+3 


593 


1 + 0.4 


Suzaku (2008/03) 


XIS 0+1+3 


592 


8 + 0.4 



Error is 1-a of the Gaussian of the folded result (e.g. figure 5b). 
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be a foreground source, not a member of the GC sources. Uchiyama, Y., et al. 2008, PASJ, 60, S35 

Conversely, Sidoli et al. (2006) suspected that SAX 
J1748.2— 2808 is located near the GC region at about 
8.5 kpc, because it has a large absorption of 1.4 x 
10 23 H cm~ 2 . With this distance, they estimated the lumi- 
nosity to be ~ 10 34 erg s _1 , significantly larger than any 
other magnetic CVs, and hence the authors declined to 
suggest a HMXB origin. We, however, note that a large 
fraction of the absorption of IPs is due to the circum- 
stellar gas, and the absorption value could be up to sev- 
eral xlO 23 H cm~ 2 (Ezuka & Ishida 1999). Therefore, 
the large absorption of SAX J1748. 2-2808 would also be 
due to circum-stellar gas, rather than interstellar gas in- 
tegrated along the long line of sight. This large amount of 
circum-stellar gas can be naturally explained as being the 
origin of the strong 6.4 keV line from neutral irons (see 
e.g. Ezuka & Ishida 1999). 
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